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Abstract 

We present a formulation of non-Abelian gauge theories in general axial gauges 
using a Wilsonian (or 'Exact') Renormalisation Group. No 'spurious' propagator 
divergencies are encountered in contrast to standard perturbation theory. Modified 
Ward identities, compatible with the flow equation, ensure gauge invariance of phys- 
ical Green functions. The axial gauge nA = is shown to be a fixed point under 
the flow equation. Possible non-perturbative approximation schemes and further 
applications are outlined. 
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1. General axial gauges IjJ-J/]] have enjoyed considerable attention amongst the non- 
covariant gauges, especially for computations in QCD at vanishing [0j>P| or non- vanishing 
temperature The main reason for their popularity stems from the fact that the 
ghost sector decouples. The number of Feynman diagrams in a perturbative loop expan- 
sion is reduced, leading to an important simplification from a technical point of view. 
Furthermore the problem of possible Gribov copies ||, generically present in covariant 
gauges, is absent 0. The price to pay is that the (perturbative) propagator receives 
'spurious' poles, which have to be dealt with separately. The question about how to reg- 
ularise the propagator as to allow for a consistent loop expansion stimulated extensive 
investigations 0. The intricacies concerning these regularisations partly spoil the advan- 
tage of having fewer diagrams to calculate. Nevertheless it has been an appealing gauge 
to e.g. calculate expectation values of Wilson loops which serve as order parameters for 
confinement. In the strong coupling limit they are expected to fulfil Wilson's area law 
which is correlated to a linear quark potential. The proper calculation of these expecta- 
tion values may also necessitate the inclusion of topologically non-trivial configurations 
like instantons [f|. Wilson loop calculations have also been used as a testing ground for 
the consistency of calculations in general axial gauges. 

Thus it would be interesting to see whether an alternative approach (other than stan- 
dard perturbation theory) to gauge theories in general axial gauges is available which 
preserves the above mentioned benefits without encountering 'spurious' divergencies. 

In this Letter, we shall argue that these 'spurious' poles are indeed an artifact of a 
perturbative loop expansion. The remedy we propose is known as the 'Exact' (or Wilso- 
nian) Renormalisation Group flOHFll- This approach has already been applied to scalar 
||13|| , Abelian Higgs theories [[TJ[]-[[16] and non-Abelian gauge theories ||17||-||21|| in covari- 
ant gauges, and is particularly useful in cases where perturbative expansion parameters 
tend to be large. The key difference to a perturbative loop expansion consists in the fact 
that flow equations integrate-out quantum fluctuations mode by mode while integrating 
over infinitesimal momentum shells. This connects continuously the classical with the full 
quantum effective action and can be interpreted as a 'coarse-graining' of the microscopic 
field theory. The perturbative loop expansion lacks the notion of 'coarse-graining' as all 
modes will be integrated-out in one step within a given loop order. 

For the issues addressed in this Letter fermions will act only as spectators. Thus for 
the sake of simplicity we concentrate on the pure non-Abelian gauge theories. 

2. Let us shortly review the appearance of 'spurious' propagator poles related to an 
axial gauge fixing in standard perturbation theory. We will start with the Euclidean 
action for a pure non-Abelian gauge theory, given in d dimensions by 

Sa\a] = -Jd d x f; v f; u (i) 

with the field strength tensor 

F; v = d,Al - d u Al + gf a hc A\A\ ( 2 ) 
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and the covariant derivative 



Df{A) = 6 a % + gf ch A% [t\ t c ] = f a bc t a . (3) 
A general axial gauge fixing for the (fixed) Lorentz vector is given by 

= I J d d x n,A; n u A a v . (4) 

The gauge fixing parameter £ has mass dimension —2 and may also be momentum de- 
pendent. In particular, the case £ = (£p 2 = —1) is known as the axial (planar) gauge. 
The propagator P^ v related to S = Sa + is 

*fj,v — 2 ' f \2 2 2 ' \ ' 

p z [np) z p z p z np 

It displays the usual IR poles proportional to 1/p 2 . In addition, we observe additional 
divergencies for momenta orthogonal to n M . These poles appear explicitly up to second 
order in 1/np and can even be of higher order for certain np-dependent choices of £ ||. 
For the planar gauge, the spurious divergencies appear only up to first order. 

This artifact makes the application of perturbative techniques very cumbersome as an 
additional regularisation for these spurious singularities has to be introduced. 

3. The key idea of the renormalisation group a la Wilson is the step-by-step 'integrating- 
out of degrees of freedom'. This program can be achieved simply by adding a scale- 



dependent term to the action fl4|-|pi 



AkS[A] = llw)~ d A » R ^ ip)A »- (6) 

Here we have introduced the infra-red (momentum) scale k which will interpolate from 
some UV scale A to the IR limit k = 0. Eq. ([]) is quadratic in the fields and leads 
therefore to a modification of the propagator. 

The scale dependent Schwinger functional W4[J] related to S k = Sa + Sgf + AfcS* is 
given by 

expW k [J] = JvAexp(-S k [A] + J d d x A^ , (7) 
and the scale dependent effective action Y k is defined as the Legendre transform of (0) 

r k [A] = J d d x j;a; - w k [j] - A k s[A], a; = ^Hl. (8) 

For later convenience, we have subtracted A k S from the Legendre transform of W k . The 
regulator R k enjoys the following properties: 

(i) It has a non- vanishing limit for p 2 — > 0, typically R —> k 2 . This precisely ensures the 
IR finiteness of the propagator at non-vanishing k even for vanishing momentum p. 
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(ii) It vanishes in the limit k — ■» 0. In this limit, any dependence on R k drops out and 
Tfc^o reduces to the full quantum effective action T. 

(iii) For k — > oo (or k — > A with A being some UV scale much larger than the relevant 
physical scales), R k diverges like A 2 . Thus, the saddle point approximation to (|7|) 
becomes exact and T k ^\ reduces to the gauge-fixed classical action Sa + S g f. 

As a consequence, the functional T k interpolates between the gauge-fixed classical and 
the full quantum effective action. The corresponding flow equation follows from (H) as 

^^Trjfrin^r 1 ^) 0) 



m / 

where the trace sums over all momenta and indices, t = In k, and 

p(2)o6/ a _ 5 2 Tk , in s 
k >» AX,X) ~ 6A°(x)6Al{x>Y 

Note that for any given scale k, the main contributions to the running of T k in (|3|) come 
from momenta around p 2 ~ k 2 . This is so because dtRk is peaked around p 2 ~ k 2 , and 
(exponentially) suppressed elsewhere. The physics behind this is that a change of T k due 
to a further coarse graining (i.e. the integrating-out of a thin momentum shell around k) is 
dominated by the fluctuations with momenta around k. Contributions from fluctuations 
with momenta much smaller /larger than k should be negligible. 

4. What have we gained with the propagator P k ^ v related to Sfc? Let us specify the 
regulator as 

K»M = 6ab [ r (y)p\» ~ f(y)p,Pu] (ii) 

and y = p 2 /k 2 . A typical class of regulator functions with the above properties is given 
by (m > 1) 

r m (y) = — !- . (12) 

(The limit m —> oo corresponds to the sharp cut-off limit JRI.) The propagator takes the 
form 

"k,iiv = cli— + a 2 — i — r a 3 — — r ho 4 — , [L6) 

pz p<t p z [np) n z p z 

with the dimensionless coefficients 

ai = 1/(1 + r), a 2 = (1 + r)(l + &> 2 (1 + r))/z 
a 3 = —(1 + r)s 2 /z, a 4 = — (r — r)/^ 

and 



(14) 



s 2 = (npf/{n 2 p 2 ) (15) 
z = (1 + r )[(l + f)s 2 + (r - f)(l + p 2 ^(l + r))]. (16) 
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The most convenient choice for practical purposes is r = 0. The cut-off function r obeys 
the limits (y = p 2 /k 2 ) 

limr m = y~ m , lim r m = 0. (17) 
For f = the propagator Pfc )At „ has the limits 

I / TI TI i \ 

limP^ = P M „, HmP feiP = - + j <U, (ig) 

lim P fciAti , = , lim P fc)AW , = P^, 

k 2 — >oo p 2 — >oo 

with P M1 , given by @ . The infrared regulator does not contribute for both large momenta 
or k — > 0. For k — > oo the propagator vanishes as all quantum modes are suppressed. By 
construction, the propagator (|13D is IR finite for any k > 0. The important observation 
is now the following: In contrast to the perturbative propagator P^, the limit of P*^ 
for np — » is finite! This holds true even for an arbitrary choice of £ (p, n) and leads to 

1 <W , l + f p^pv 1 n^riv 



1 + r p 2 (l + r)(r — f) (1 + r)(l + p 2 £(l + r)) n 2 p 2 



(19) 



Thus (^) is perfectly well-behaved and finite for all momenta p (as long as the regulators 
r and f have not been chosen to be identical). It is noteworthy that the 'spurious' 
divergencies are absent as soon as the infra-red behaviour of the propagator is under 
control. 

5. Now we shall discuss the implications of gauge invariance, and in particular the 
question about how to ensure gauge invariance of physical Green functions. To that end 
we consider the Schwinger functional as given in ([?]). Invariance of the measure DA 
under the infinitesimal gauge transformation A — > A + Da leads to the so-called modified 
Ward-identity (mWI) 

= (Dpi - n,Df{x) -L n.A^x) - Df{x) R^ v A c u (x)^ . (20) 

The mWI contains quadratic powers of the gauge field, which is a consequence of the 
gauge fixing and the cut-off term. Usually one converts the WI into a linear form by 
introducing ghosts and using BRST-invariance. However, due to the cut-off term this is 
no longer possible. With 

atpM = OfS? + Rt,y l (x,x') (21) 

the mWI can be converted into 



„, ( A STk[A] n „,, ; . 1 ,,,, 

^ ■ (22) 



WftA] = D;\x)-^-± - n,Df{x) -^n„Al(x 



5A b ^{x) ^ n 2 C, 

J dd V f abc i^i n ^ 6Cd + R tr) (*,y)Gt,(y,x) 
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The first term is the usual covariant functional derivative of T^, the second one stems 
from the gauge fixing, while the third one, a consequence of the regularisation and a 
possible momentum dependence of £, introduces non-local contributions. Note that the 
.Rfc-dependent term vanishes for both R k — > oo and R k — > 0. Both ( [201) and (|2~2"D reduce to 
the usual Ward identity in the limit k — > [0. For k ^ 0, they explicitly depend on the 
regulator, and the relations amongst different n-point functions of the theory as implied 
by ( p2|) will depend on the particular form of the coarse-graining. 

We will now derive the important result which is that the flow equation flj^) and the 
mWI ( p2]) are compatible. To that end we calculate the scale dependence of (0). Using 
the flow equation (Q) for r&, one can check explicitly that 

= > ( G *^n x w *° (23) 

where the trace sums over momenta and internal indices. It follows that if I\ fulfils the 
Ward identity at some scale k Q , and if it evolves according to the flow equation (|9|), then 
the new action will again fulfil the corresponding mWI.Q Thus it is in principle sufficient 
to fulfil the mWI at some initial scale A in order to ensure gauge invariance of the physical 
Green functions, that is in the limit k — > 0. 

Let us comment on the possible momentum dependence of £. As £ is dimensionful, it 
is natural to consider it as an operator. In perturbation theory, for example, the planar 
gauge £(p) = — 1/p 2 introduces a momentum dependence in order to reduce the degree of 
the spurious divergencies. This is no longer necessary as the divergencies are absent in 
the present approach. Therefore we may restrict ourselves to a momentum-independent 



gauge parameter. This is an important simplification, as the mWI (22) reduces to 



D ^ {x) li§) ~ nvK{x) - g j ddy f abCR t^ y) G t»(y> x ) = °- (24) 

As a consequence the possible tensor structure of vertices is considerably simplified. How- 
ever, one may ask whether the limit £ — > defines the axial gauge also for momentum 
dependent £(p). In perturbation theory, the answer is yes p. In the present context, it 
suffices to show that the term 

j^y^^ft,M (25) 

in ( p2"D vanishes in the limit £(p) — > 0. That this is indeed the case can be seen as follows: 
Expanding in powers of the gauge field yields expressions containing the effective 
propagator and n-point vertices T^ >2 \ Both G° k and r^ 2 -* are completely determined 
via the (finite) flow equation. Thus they cannot exhibit singular terms proportional to 

3 It can be shown that ( p3| ) is valid for general linear gauges. For covariant linear gauges, a similar 
result was obtained in fllSl. 
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l/£. Using in addition n^Gk^u = 0(0 (due to the gauge fixing) we find that ( p5| ) is at 
least of order 0(£) for £ — > 0, which establishes the axial gauge even for a momentum 
dependent gauge fixing. 

6. As a first application of this formalism we shall argue that the choice £ = corre- 
sponds to a fixed point w.r.t. the flow equation. In principle, this can be shown through 
an explicit computation of the flow equation for £(&). However, this would necessitate 
an explicit Ansatz for (and would therefore be only an approximation). In contrast, 
we present an argument which makes use only of the mWI. With the preceding result 
in mind we can safely restrict ourselves to the case of £ being momentum independent, 
although the result will hold true for general £(p). First note that £ enters the mWI both 
explicitly and implicitly. The £ appearing explicitly corresponds to the choice of £ at 
some initial scale A, £ = £(A). An implicit dependence occurs through the dependence 
of T fc on the scale dependent £(&). Let us choose £(A) = with I\ solving ( |2"2"| ) and 
assume that £(&) ^ for some k < A. This means in particular that will no longer 
contain a singular term ~ l/£. Thus the only singular term appearing in the mWI is the 
term explicitly proportional to l/£. One can always find an A a such that nd nA a does 
not vanish. Therefore with £(&) ^ can not be a solution of (|22|) for £(A) = 0. But 
this cannot be true as the compatibility of the flow equation and the mWI fl23|) implies 
that Tfc solves the mWI. It follows that £(&) = for £(A) = 0. Hence the axial gauge is 
indeed a fixed point of the flow equation. Note that this argument does not involve any 
approximations regarding IV 

We conclude that out of all general axial gauges (even momentum dependent ones) 
the axial (£ = 0) gauge is singled out and appears to be the natural choice. Furthermore, 
£ = is well-suited for actual computations as both the flow equation and the mWI are 
rather simple in that case. 

7. A second application concerns the question about how to control gauge invariance 
for an approximate solution of the flow equation. Generally speaking, the task of solving 
the flow equation for gauge theories faces two main problems. The first one is to find 
a solution of (|22| ) at some initial scale k = A in order to ensure gauge invariance of the 
physical Green function in the limit k — > 0. The second one concerns an Ansatz for the 
functional form of T^A]. As not all possible operators can be taken into account, the 
validity of (]22[) for k < A (automatically ensured only for the full effective action T^A]) 
can no longer be taken for granted and has to be checked independently. 

As an example, consider the action r = Sa + S g t for £ = at the scale k = 0. Obvi- 
ously, To is a solution of the mWI fl24]) . The flow equation can be integrated analytically 
in leading order of perturbation theory, replacing T k through T on the r.h.s. of (0), to 
give 

r fc = T + ^Tr In (Tf + R k ) - In T^. (26) 
The mWI without the ^-dependent term is solved by r . For the remaining terms in 
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(M) we obtain in leading order 



±D*{x) S{Tr ln(r » 2) +*f x ~ ^ ln F ° 2)) - g j d d y r b ^(x, y)G* „(y, x) = 0. (27) 



Tq is the second derivative of T with respect to the gauge field fllPl) . We use the 

<5A' (SA) 2 



commutator [Djj, tt^W] and (]23) for T to obtain 



^'' ""S^ = - g I dS r*^^ Urirf m - (28) 



For k — 0, (p8|) is simply zero. Using (To 2 " 1 + i?fc) 1 = Gfc + 0(g) and inserting (|28|) into 
( p7|) results in 

-0 / A r" c (rg? + (x, y)Gj, M (y, x) ~ -g 2 f abc 6 bc = 0. (29) 

Thus we have shown that the compatibility of the flow and the mWI can be maintained 
even within an approximate solution. It is straightforward to show that this holds true 



systematically even for higher orders within a perturbative loop expansion |21| . 

Note that in jTjl] the 1-loop perturbative compatibility of the flow equation was checked 
explicitly (for covariant gauges) for the scale dependent gluon mass parameter using a 
BRST-formulation. This is a rather non-trivial task since the flow equation as computed 
directly from ([|) or from the corresponding Slavnov- Taylor identity receives contribu- 
tions from quite different diagrams (involving ghosts and gauge fields). However, in our 
formulation (without ghosts) the consistency check is rather simple and is done without 
problems for the entire effective action. 

In more general situations, and especially within non-perturbative regions, it is not 
obvious how the compatibility between flow and mWI of a given truncation can be main- 
tained. However, it is still possible to exploit the compatibility condition and to use it 
as a control mechanism for the Ansatz itself. A natural implementation would be to use 
the mWI as a flow equation for some of the relevant couplings (like mass terms) of the 
theory. Comparing the flow of these operators with the flow as derived directly from (|9]) 
allows one to control the domain of validity of a given truncation.^ 

For approximations beyond perturbation theory the flow equation @ and the mWI 
(p4]) can also be used to control the dependence on of the effective action and thus 
generalise the observations of to the case with a cut-off term. 

8. Let us finally comment on the computation of the 1-loop /3-function.Q This is 
another crucial test for the viability of this method. In standard perturbation theory, 
the 'spurious' singularities seem to play an important role and do already contribute on 



4 See |lq| , [19| where a similar line of reasoning has been employed on QCD in covariant gauges. 
5 All the details of the computation shall be presented elsewhere p3j- 
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the 1-loop level 0. In the present context they are absent throughout and it would be 
interesting to check that the correct 1-loop running is still coming out. 

The running gauge coupling g\ is related to the scale dependent wave function renor- 
malisation of the field strength Zp k via g\ = k^^g 2 / Zp^. Using the Ansatz 

r k = \J d d x z F , k f; u f; u + s gi [z l J 2 k A] (so) 

and (H) one can deduce the 1-loop flow for Z F ^ k by projecting on the F 2 -terms in the flow 
equation. This is self-consistent at the 1-loop level. The traces involved in (|9]) can be 
calculated using heat kernel methods^ and one obtains 

dt \nZ F , k = ^g 2 + 0(g'). (31) 



The /9-function follows immediately from (|3l| ) as 

fa=d t gl = -^gi + 0(gl), (32) 

i.e. the well-known 1-loop result. The above (universal) result can be shown to hold for 
any regulator with the properties (i) - (iii). 

9. To sum up, we have developed a self-contained and self-consistent formulation of 
QCD in general axial gauges which allows full control over the IR behaviour of the theory. 
In contrast to standard perturbation theory the 'spurious' propagator singularities are 
naturally absent. Gauge invariance of physical Green functions is controlled via the mWI 
which is shown to be compatible with the flow equation. The absence of ghosts results in a 
rather simple expression for the m WI. We have shown that of all general axial gauges the 
axial gauge nA = is singled out as it corresponds to a fixed point under the Wilsonian 
flow. Possible ways of finding approximate solutions even beyond perturbation theory 
have been briefly indicated. The calculation of the 1-loop running of the gauge coupling 
has been outlined and the /3-function agrees with the known perturbative result. 

This formulation seems to be a good starting point for theories where Lorentz covari- 
ance is naturally broken. This is the case for QCD at finite temperature where the heat 
bath singles out a rest frame. It is straightforward to apply the present approach within 
the imaginary time formalism. The propagator would still come out without any 'spuri- 
ous' divergencies, in contrast to the recent proposal based on a renormalisation group 
for the temperature fluctuations only. Another possible application concerns quantum 
field theories with both electrically and magnetically charged U(l)-fields, as this necessi- 
tates the introduction of a fixed Lorentz vector [^|. In either case, the Lorentz vector 
should be used for the axial gauge fixing as described above. It may also necessitate the 
inclusion of topologically non-trivial configurations into the effective average action |20[ . 
We hope to report on these matters in the future. 

6 Here, the heat kernel is not used as a regularisation since (M) is finite anyhow. 
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